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In this paper, a variational perturbation scheme for nonrelativistic many-Fermion systems is 
generalized to a Bosonic system. By calculating the free energy of an anharmonic oscillator model, 
we investigated this variational expansion scheme for its efficiency. Using the modified Feynman 
rules for the diagrams, we obtained the analytical expression of the free energy up to the fourth 
order. Our numerical results at various orders are compared with the exact and other relevant 
results. 



I. INTRODUCTION 



For calculating the free energy of a system, there exist usually two basic methods, the conventional perturbative 

and the variational methods [n|. In fact, these two methods are standard in calculations of many physics problems. 

m However, it is well known that the former is useful only for small perturbing potentials, whereas the latter lacks 

^ systematic schemes to control its accuracy, albeit it is valid for any potential. In order to overcome these difficulties 

*■ and improve the variational method, a variational perturbation idea of properly combining the two methods was 

^-v pioneered by Koehler in lattice dynamics and Seznec and Zinn- Justin |^ on an anharmonic osillator decades ago. 

^^ Later, the idea was extended to path integrals by Feynman and Kleinert Q as well as Okopinska Q], and further 

f^ developed by other authors Mm Q Very recently, three of the authors (You, Kim and Nahm) and their collaborator 

^^ presented a variational perturbation scheme for a nonrelativistic many-Fermion system in the functional integral 

^^ formalism H . In order to test the efficiency of the method, Ref. H provided a model numerical calculation at the 

^ zero temperature up to the second order. Obviously, this test is limited, and a more wide investigation is necessary 

O f or the variational perturbation scheme. 

. K . Actually, the scheme in Ref. |q] is a Taylor series expansion based on the variational result obtained in the spirit 
^ of Feynman variational principle fit]. This scheme can improve the variational method if used judicially. We note 
Cu that there exists no work performing the same scheme for Bosonic systems. Although Okopinska H] and Krzyweck 
p~ M established two kinds of expansion schemes, the optimized expansions and the cumulant expansions respectively, 
.. yvhose lowest order results are the variational result, they are really not based on the variational result because 
^ the variational procedure was performed at the truncated order. Further, improvement to the variational method is 
k> notoriously difficult, and different schemes will have their own advantages over others. So it is worthwhile generalizing 
''T^ the scheme in Ref. g to Bosonic systems. 

C^ ' For the above two purposes, an anharmonic oscillator will be an effective laboratory. For a one-dimensional 
anharmonic oscillator, the Hamiltonian is 

H{t)^^p^ + lmuj^x^{t) + Xx\t) . (1) 

2m 2 

Here, x is the space position, p the momentum, t the time, and m, lu as well as A are the mass, frequency and coupling 
strength, respectively. Such an anharmonic oscillator is probably the simplest model which does not have exact 
analytic solution. Its exact free energy was obtained numerically |1Q,p1. Moreover, there exist many approximation 
methods to calculate its free energy p|,|6|,p|,[ll| . All of these make Eq.(l) an ideal candidate for our aim. Therefore, in 
this paper, we will generalize the scheme in Ref. [||, taking the simple anharmonic oscillator, Eq.(l), as a laboratory 
and investigate the Taylor expansion scheme on the free energy based on the variational result obtained from the 
Feynman variational principle m . 



*permanent address, E-mail: wenfalu@online.sh.cn 
^Here we are far from exausting the relevant literature. 



Simultaneously, the anharmonic oscillator itself is useful in chemical physics ]1^ and many physics problems, such 
as thermal expansion, phonon softening and structural phase transitions [ |l3| . Although its free energy was calculated 
by using many methods, their numerical (approximate and exact) results ^ were focused mainly on the case of the 
reduced temperature T less than 1, except for Ref. [Q which provided accurate results at moderate T. Therefore, 
we also test whether our scheme can produce reliable results at moderate temperature range. 

We will work within the functional integral formalism [Q . Next section generalizes the scheme in Ref. g to calculate 
the free energy of a system with the potential V{x). In Sect. Ill, we carry out the scheme on the anharmonic oscillator, 
Eq.(l). Sect. IV presents calculations of the free energy for Eq.(l) up to the fourth order and makes comparisons with 
the exact, and various approximate results, such as, variational, cumulant-expansion and optimized-expansion results. 
Conclusions are given in Sect.V. 

II. A VARIATIONAL EXPASION OF FREE ENERGIES FOR BOSONIC SYSTEMS 

For a system with the Hamiltonian H(t) = -^p^ + V"[a;(t)] with y[x(i)] = ^muj'^x'^{t) + Vi[x{t)], the generating 
functional is Q (1990) 

Z[J] = f V[x{t)] exp{- / [^mx{T){-d',)x{T) + V[x{t)] - J{T)x{r)]dT} , (2) 

Jx{a)=x{l3} Jo ^ 

where t = it is the imaginary time, 9,- = ^, /3 = -^ with k the Boltzmann constant (we will consider it as unity 
for convenience) and T the temperature. J is an arbitrary external source and the symbol I?[a;(T)] represents the 
functional measure. In Eq.(2), x{0) = x{(3) means that the functional integral should be executed over all the closed 
paths |1J] (1990). To introduce a variational parameter 17, one can rewrite the classical action functional S\x\ in the 
Euchdean space-time in Eq.(2) as |0,§,3 



'^ 1 1 



S[x, J] - / [-mx(T)(-a^ + n^)x{T) - J{t)x{t) ~ -mn^x\T) + V[x{t)]]A 



T 



^2 ^ '^ ^ ' ^ ' ^ ' ^ ' 2 

= Sa[x,J]+SD[x] (3) 

withS'o[x,J] ^ jl^[\mx{T){-dl + ^'^)x{T)~ J{T)x{T)]dT ixnd Sd[x\ = j^[-\mn'^x^{T) + V[x{T)]]dT. Thus, Z[J] can 
be rewritten as 



Z[J] = [ V[x{t)] exp{~So[x, J] - Sd[x]} 

Jx(0)=x{P) 

= exp{- / [-\mnH\ + V[5J^]dT} [ :D[x(r)] exp{~So[x, J]} 

Jo ^ Jx{0)=x{P} 

r^ 1 1 

= (i?ei(G-i))-5exp{- / [--min" -u:-')5]^+Vi[5j^]dT}exp{-JrGrr'Jr'dTdT'} , (4) 

Jo ^ ^ 

where, Jr = JIt), Sj^ = gj-, G^^ represents the operator —9^ + fl^ with the propagator Gtt' — G{t,t') and Det 
means the determinant. In the third equation of Eq.(4), we have carried out the Gaussian functional integration Mm 
(1990). Note that Grr' can be expanded into a series owing to the closed path requirement in the functional integral, 
Eq.(2) |14| , |l5| . To calculate the free energy, we express the partition function Z = Z[J ~ 0] in the following form, 

^f ^r ^_M„ r cr , _ nn ^-(o)=-(/3) ^[^^^^1 '"''P^^^ot''' '^ ^ °] " '^^t'']^ 



r 

/ I?fx(T)l eXpj — 5o[x, J = oil- „ „r / M r ^ r T „n 

Jxio)=xiP) ^ ^ '^ ^ ^ ^^ X.(o)^^(^)P[x(r)]exp{-5o[a;,J = 0]} 

- {Det{G-^))-^ < exp{-y [-^m{n^ - lj^)SI +Vi[Sj^]]dT} >g 

r^ 1 

= (I?e<(G-i))-5 exp{- / < --m(f72 _ ^2)^2^ ^ y^^^^^^ >^ ^^| 
Jo ^ 



(5) 



In this paper, an anharmonic oscillator does not include the double-well potential case. 



< exp{- ^ [-^rn(f72 - uj^)6% + F,[<5jJ 



< -im(r!2 - uj^)5\ + ^/[(JjJ >G]dT} >G . (6) 



Here, we have used the fohowing notation and relation 



U)=.IP) n^ir)]0[x] exp{- J^ \mx{T){-dl + n^)x{T)dT} 
< 0[x\ >G ^ 



=< 0[6.j] >G . (7) 



— 0[Sj]e:xjp{-JrGrT'JT'dTdT'} 



J=Q 



Obviously, when V/[a;(t)] is not zero, it will be impossible to obtain analytically exact partition function and, hence, 
one has to design some scheme to produce an approximate solution. For the case of small Vj, one can take il = lu, 
make a Taylor series expansion to the exponential in Eq.(5) and then truncate the series at some order to approximate 
Z . This is just the conventional perturbation theory and the propergator Gtt' is the bare propagator. When Vj is not 
so small, the above perturbation method is valid no longer. In such a case, an effective alternative is the variational 
method which is based on the Fenyman variational principle. In the following, we briefly introduce the variational 
method. 

Exploiting the Jensen's inequality ||l|,^, one can have 

< exp{- J [-^m{n^ - iu^)6l + Vj[Sj^]]dT} >g 

> exp{- / < --m{n^ - LU^)5X + Vi[dj^] >g dr} . (8) 

Jo ^ 

Substituting the above equation into Eq.(6) leads to a relation for the lower limit of the partition function, i.e., 

/■^ 1 

Z > {Det{G-^))-^ exp{- / < --fn{n^ - uj^)5j^ + Vi[6j^] >g dr} . (9) 

Jo 2 

Hence, the free energy is 

F^-^HZ)<^HDet{G-')) + ^J <-lrnin'~Lo')Sl+Vi[Sj^]>GdT^F. (10) 

Obviously, making F of the last equation the absolute minimum will lead to a minimum upper limit of the free energy, 
Fq, with the variationally extremized condition, 

§,-0. (11) 



and the stabilized condition 

S^F 



{sn^y 



>0. (12) 



The parameter fi which renders F absolutely minimized will be chosen from the three possibilities; the non-zero 
solution of Eq.(ll), zero and oo. F with such an Q is just Fq, the variational result of F. This procedure is essentially 
same as was done in Ref. [|l6| for Eq.(l). 

Entering the above variational result into Eq.(6) and taking the logarithm, we obtain the following expression of F: 

F^Fo -Ih< eM-[SD[Sj]- < Sd[Sj] >g]} >g] ■ (13) 

P 

Now, we make a Taylor series expansion of the exponential in the last equation, and the average < • • ■ >g can be 
calculated order by order through borrowing the Feynman diagram technique [|l4| Q The logarithmic operation in 
Eq.(13) is equivalent to discarding disconnected diagrams [Q. Consequently, we have 



^Now diagrams are no longer the bare Feynman diagrams in the sense of perturbation theory owing to Q, y^ cu which has been 
variationally determined. 



F = Fo+^F(") (14) 

2 

with the nth order correction to the variational result 

p(n) ^ (-l)"+lil < [Sn[Sj]\r^r,- < SD[5j]\r=r, >g] 

■ ■ ■[SoiS.Alr^r.- < SD[Sj]\r=r. >g] 

■ ■ ■ [^D[^j]|r=r„- < SD[Sj]\r=r,, >g] >G,C ■ (15) 

Here, the subscript C means that only the connected diagrams have their contributions to the free energy ||lj]. Eq.(14) 
corresponds to a systematic Feynman-diagram-like expansion, and from it, one can estimate the approximate values 
of the free energy, F, order by order. Thus, we have finished the generalization of the scheme in Rcf. [g[ to a Bosonic 
case. In this scheme, the parameter f2 is variationally determined before the series expansion is performed, and it is 
identical for all orders. Obviously, this is a Taylor series expansion around the variational result, and so we call it the 
variational expansion. Because the ith factor in Eq.(15) (i = 1,2, • • -n) has the term — < SD[Sj]\T=Ti >g= ~ /q < 
— ^rnil^S^ + V[Sj^, ] >G dTi, which has negative sign against the major part of Fq, one can expect to get a simplified 
diagram rule, as will be shown for the system, Eq.(l), in the next section. Next, we apply the above procedure to the 
anharmonic oscillator and make a comparison with existing results in the literature. 

III. APPLICATION TO THE ANHARMONIC OSCILLATOR 

For the system, Eq.(l) (F[x(i)] = ^■muj'^x^{t) + Xx^{t)), the procedure from Eq.(8) to Eq.(12) yields easily the 
variational free energy, Fq, 



with the variationally extrcmized condition {j^ = 0) 

^'--' + ^cotH^). (17) 

Here, owing to the periodicity of the path in Eq.(2), we have used the following propagator 

^- pl-^rniu^^ + n2f 2mn sinh(^) ^ ' 

with ujn the Matsubara frequency |^,|5[. Eq.(16) coupled with Eq.(17) are just the variational result of F in Rcf. |16[ |. 
Using the relation ^m'^{uj'^ ~ il^) = — GAGrr from Eq.(17), we have 

< Sd[6j] >g = < ^Mij^ ~ ^^)SX + XSX >G dT 
Jo ^ 
fP I rP 

= / [-to(cj2 - n^)Grr + "iXGrrGrrW = SXGrrGrrdr . (19) 

Jo 2 Jo 

So, for any i, one has 

< • • ■ SD[Sj]\r=T,'- < Sd[Sj]\t=t, >g] ■ ■ ■ >G,G 

r^ 1 1 

= < • • • y [-to(l^2 _ ^2)^2^^ ^ ^^4^^ _ ^ --m{uj^ - n')6\ + \5^j^^ >G dn--- >g,g 

rP 1 
= <■■■] [^rniuj-" - n'){{5.,^f + G,,,J 

+ X{C5.uX + QGr.rXSj.f + iGr,r,Gr,r,) + 3AG,, ,, G.^.Jdr, • • ■ >G,G 

= <■■■ f\x{C6jJ*--->G,G, (20) 







...... 1 

where, the symbol o in J means that the functional derivative with the index i takes effects on expl-^JrGrT'Jr'dTdT'} 
only if it makes up pair with any other functional derivative with the index j ^ i to yield Gnr ■ In going to the second 
step of the last equation, a concrete analysis has led to the following equivalent properties (o means equivalence): 

S^ ^ Csjr, )' + G.,.. , Sj^^ ^ Csj.^ )' + (^Gr^r, Csj., f + 3a,nG,,r. . (21) 

In terms of the Feynman diagram language, it implies that only the legs which come from different vertices can connect 
each other. 

Substituting Eq.(20) into Eq.(14), we can estimate the higher-order corrections to Fq in Eq.(16) with the help of 
the diagram technique |14| . The free energy for Eq.(l) is now 



P Jo n 



(22) 



with the nth order correction, 



1 A" 



/9 „ rf^ 



P "! Jo Jq 

Here, the modified Feynman's rules for drawing diagrams are quite simple, and they are as follows: 

(1). . . Propagator, G^ra ; (2). \^ Vertex, '^^Jq dr . 




For the nth order, there is an additional total factor —75^ ■ From Eq.(23), it is evident that there will be not any 
Cactus diagrams appearing at any higher order, which is demonstrated by the diagrams in the next section. This 
simplifying feature of diagrams is similar to what occurs in the Fermionic case [pi . A further analysis indicates that 
there exist the following four types of building bricks for any nth-order connected diagrams (n > 2) : 







(a) (b) (c) (d) 

which correspond to the four kinds of partitionings of the integer "4" : (a) 2-f 2 , (b) 2+1+1 , (c) 3+1 and (d) 
1+1+1+1, respectively. In this figure, the intermediate vertex of the brick (a) has two legs connected with one 
different vertex and the other two legs with another different vertex; the left (or right) vertex of the brick (b) has two 
legs connected with one different vertex and the other two legs of it will be connected with some other two different 
vertices respectively; the left (or right) vertex of the brick (c) has three legs connected with one different vertex and 
the other leg of it will be connected with some different vertex; and the vertex of the brick (d) will have its legs 
connected with some four different vertices respectively. These four bricks are helpful for drawing various distinct 
diagrams at any order as one can see from the five diagrams drawn in the next section. For example, all of them do 
not contain the brick (d), the second-order diagram consists of only the brick (a), so does the first diagram of the 
forth-order diagrams (4a). 

In the next section, we calculate the free energy up to the fourth order from Eq.(22). 

IV. ANALYTICAL EXPRESSIONS AND NUMERICAL RESULTS 
UP TO THE FOURTH ORDER 

According to the last section, the topologically non-equivalent diagrams at the second, third and fourth orders can 
be drawn as follows. 







and 




2nd 



3rd 



4th (4a) 



4th (4b) 



4th (4c) 



For these five 
diagrams from the left to the right, their symmetry factors (the number of topologically-equivalent diagrams appearing 
in the expansion) are N2 = 4!, iVg = |^ • (2 • Clf , A^4a = :^ • (2 • Clf, iV4h = ^ • (C| • 2 • Clf ■ 2^ and 
7V4c = ^ • (C4 • 3! • C^Y ' 2, respectively. Thus, one can easily write down the corrections F^^^ F*^"^' and F^''^ 
according to the above diagrams and then calculate them as 



pm 



^^sinh-4(^)[6/3f7 + 8sinh(/3f7) 



sinh(2/3f7)] 



(24) 



i^(3) 



1A3 



-N. 



dT.dr.drsGi^^^Gi^^Gi^^,^ 



/3 3! 

+48 cosh(2/3J7) + 3 cosh(3/3f}) + 108/317 sinh(/3r2)} , 



3 + Sp^n^] cosh(/3rj) 



(25) 



and 



Jo 



/3 4! 

3A4 



2 

T4T1 



^2 



-rir2^r3r4^'^2-r3*^r2r4*^'^i'^3^Tir4 "^ -^^4c^rir2 ^r3r4*^'^2-r3 '-^T4Ti 



'3 /^3 /^ /^ 1 

'-J-r-^-i-. ^ToT-:* '"'"T/ITi J 

sinh"*(^){6291 ~ 181320/3^172 + 25920/3''r2'' 



32768^77181712 ^ 2 

+6[71 + 13156/3^172 + 2688^^17^] cosh(/317) 
+48[-134 + 2115/3^172 + 6/3^17''] cosh(2/317) 

-432 cosh(3/317) + 864/3^17^ cosh(3/317) + 141 cosh(4/317) + 6 cosh(5/317) 
-191394/317 sinh(/317) + 129456/3^173 sinh(/317) + 42568/317 sinh(2/317) 
+12672/33l73sinh(2/317) +37750/317 sinh(3/317) + 1600/317 sinh(4/317)} . 



(26) 



These analytical expressions of Eqs.(24), (25) and (26) are the main results in this section. In order to obtain them, 
we have had to handle the absolute value symbol in the expression of Gt-t' (see Eq.(18)). It is straightforward to 
calculate the integrals in F^^) a^fj ^(3) j-jy dividing the integration domains into 2! and 3! parts respectively. As for 
F^'^>(n > 4), multi-dimensional integration domain, which exceeds our direct intuition, is involved. However, for any 
77-dimensional integration domain, one can divide it into 77! sub-domains so that, for each sub-domain, the relation 
Til < Ti2 < Tij < • • • < Tj. • • • < Ti^g < Ti^_2 < ''■j„_i < Ti^^ holds. Thcu, mimicking the calculation of F^^\ one can 
find the following equivalent relation 



dri 



^ 



E 



dn. 



dr,;, 



dTi, 



dn. 



dn. 



dr. 



(27) 



where the letter "P" below the summation symbol means the summation is carried out over all the 77! sub-domains. 
Eq.(27) allows one to obtain Eq.(26) with the aid of the computer software Mathematica. 

Using the above results, we can now readily calculate the free energy up to the fourth order : F2 = Fq + F^^) , 
F3 = Fq + F(2) + F^^' and F4 = Fq + F^^) + F^^'> + F^''' . In the following, we will numerically compare them with 
existing results to examine rediability of our scheme. 

First, we compare our results with the exact results obtained from Ref. 0. Using F = ^2^1n[X]n e^'^"'^'^] (En is 
the ?7th eigenenergy of Eq.(l)), letting 777 = cj = 1, and for F < 1, one can calculate the exact free energies from 



Table V in Ref. Q. For this case, we plot Fig.l with A = 1. In Fig.l, the dotted, short-dashed, medium-dashed, 
long-dashed and solid curves are the exact free energy Fexa^ Fq, F2, F3 and F4, respectively. Fig.l indicates that : (i) 
when the temperature is near zero, F2 and F3 are very to close Fg^a, whereas F4 is unbounded from below; (ii) when 
the temperature is greater than 0.5 or so, F2 and F4 provide substantial corrections to Fq, and F4 gives better than F2 
does, while F3 is close to Fq. Here, we note that the invalidity of F4 at very low temperature is not unexpected. Since 
the present scheme is basically the Taylor expansion of the free energy, the smallness of the temperature prevents and 
competes with the convergence process of the perturbation and finally wins over at the fourth order. 

Then, we can compare our results with the accurate free energies, Faccu from the Okopinska's optimized variational 
method pM . In order to compare with Okopinska's data, we used the definitions of the dimensionless quantities in 
Ref. @, that is, m = 1, z = \oj'^\~i, ^\~^ -^ ^, TX'^ -> T and FiX'^ -> F, |. For the case of z = 10, which 
corresponds to A = 0.01118 in the dimensionlized system of Ref. |g[, we give the comparison in Table I Q From this 
table, one can see that F4 has a better agreement with Faccu than F2, F^ and Fq except for T = 30. 

Thirdly, to show the improvement of -Fb by higher-order corrections, for the range of 1 < T < 50, we plotted the 
results in Fig. 2 to compare i^4 with F^ in the cases oi z — 0.2, 1, 10,30 and 50. In Fig. 2, we use the same curve type 
to represent F4 and Fg, and between the curves within the same type, F4 is always the lower. Also, Fig. 2 shows that 
F4 coincides almost with i^o for z = 30 and 50, and the differences between Fq and F4 are quite large for both z = 1 
and z = 0.2. From Fig. 2, we learn that : (i) for a given temperature, with the increase of z, i.e., with weakening 
the coupling, the corrections of F4 to Fq get smaller; (ii) for a given z, with the increase of T, the corrections of 
F/^ to -Fb become larger; (iii) the quite large differences between Fq and F4 imply that our scheme becomes invalid 
with decreasing z or with increasing coupling strength (for a fixed w^). The third point is similar to the optimized 
expansion g (1994). 

Finally, taking m = u ~ 1, we compare our results with those obtained from the cumulant expansions [|6| in Table 
II. In Table II, the free energies F2 and F3 are ours, the free energies F/^'' and F^"^ are the first- and third-order 
results from the cumulant expansion in Ref. |^ (i.e., Fi and F2 of Table II in Ref. |^). The temperatures in this table 
is lower than 1 and the corresponding z is small. So we did not include F4 in the table due to its invalidity. This table 
indicates that our F2 and F3 are nearer to the exact value than F-^^ , but is not so good as F^"^ . This reflects that the 
convergency of our expansion is not so fast as in the cumulant expansions in Ref. 0. As for the optimized expansions, 
Ref. compared the free energy with the exact results in the case of both w = (i.e., the smallest z) and the reduced 
temperature less than 1, and demonstrated a fast convergence. Here, in Fig. 3, our results are compared at various 
orders with Fig. 2 in Ref. [pi (1994) n In Fig. 3, from the upper to lower (for large values of /?), the first, second, sixth 
and seventh curves are the variational, the third-, the second- and the fourth-order results of our scheme, respectively. 
The third and the fourth curves are the second- and the third-order results of the optimized expansions, and the fifth 
curve is the exact result. The variational result is just the first-order result in the optimized expansions. From this 
figure, our third-order result is not so good as the results in the optimized expansions, and our second-order result 
is almost as good as the second-order results in the optimized expansions. This figure simultaniously indicates the 
invalidity of our fourth-order result. Thus, the optimized expansion has also a better convergency than our results. 
Additionally, the optimized expansions approach the exact result monotonously in orders, whereas our results oscillate 
with orders. However, we want to emphasize that our sheme is not so complicated as the optimized and cumulant 
expansions. The crucial difference between ours and the optimized or the cumulant expansions is that our scheme 
performs the variational procedure at the lowest order and, accordingly, the parameter 17 is identical for all orders, 
whereas in both the optimized and the cumulant expansions PJq], the variational procedure are performed at their 
truncated order, and, consequently, 17 at one order is different from the next. It is this difference that gives rise to 
the simplicity and slow convergency of our scheme and the fast convergency and complication of the optimized or 
cumulant expansions. 

V. CONCLUSION AND DISCUSSION 

In this paper, we have generalized the scheme in Ref. [B to a Bosonic case, and taking the anharmonic oscillator, 
Eq.(l), as a laboratory, provided a wide test of its efficiency. Our investigations demonstrate that the present scheme 
can substantially improve the variational result even in the second order and when the reduced temperature is greater 



''in Ref. ||lO|, there is a typo on the rescaling expression of T, and here it is corrected. 
^The data of the accurate free energies were provided by Okopinska, the author of Ref. 
^The data from the optimized expansions were also provided by Okopinska. 



than 0.5 or so and the reduced coupUng parameter A is not strong (or z is small), the free energy for Eq.(l) up to the 
fourth order in our expansion gives a good agreement with the accurate result. We also demonstrate that for the free 
energy of Eq.(l), when the reduced temperature approach zero, or the reduced coupling A is strong (or z is small), 
the fourth-order result is invalid. Thus, from our investigations here, one can see the efhciency and limitations of our 
scheme. Here, we also note that the present scheme is much simpler than the optimized and cumulant expansions, 
albeit it is not so fast convergent as they are. We believe that a simple scheme is often necessary and useful because 
the non-perturbative method beyond the Gaussian approximation is extremely complicated in general. Additionally, 
we gave the approximate free energy of the system, Eq.(l), for moderate temperature range. Although the exact 
results for the moderate temperature has existed in the literature lin], our results can be readily used, as a basis of 
quantitative comparison, when some other approximate methods produce the free energy at the same temperature 
range. 

In general, it should be noted that the variational perturbation theory yields an asymptotic rather than a convergent 
series H] , and hence a concrete range of its validity in a specific problem has no universality. As for any specific 
problem, the present scheme should always be used with a judicial examination of the true physical property. We 
believe that the present paper can provide a qualitative reference for an application of our scheme. In particular, 
when a specific problem is too complicated to treat beyond the Gaussian approximation with other expansions, we 
believe that our scheme can provide a simple and viable tool to treat it. 

Finally, we want to point out that, although we treated the quantum-mechanical anharmonic oscillator only in this 
paper, it is straightforward to apply our method to finite temperture scalar field theory |15| . Especially, when it is 
generalized to the 0^ models |0, we expect that the simplicity of the method still holds there. 
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Fig.l For the case of T < 1 and A = 1, F2,Fs and F4 the free energies up to the second, the third and the 
fourth orders are compared with the variational result Fq and the exact free energy. The exact results were calculated 
according to F = ^Tlii[J2n e^^"^'^] and Table V in Ref. Q. Here, En represents the nth eigenenergy for the system 
Eq.(l). We took m = uj = I. In this figure, when T > 0.6, the curve for F3 almost coincides with the curve for Fq. 

Fig. 2 For the case of T > 1 and for several values of z, F4 are compared with Fq. We use the same type of curves 
to represent F4 and Fq, and the latter is always above the former. But for the cases of z = 30 and 50, F4 almost 
coincides with Fq, and for the case of z = 1.0 and 0.2, F4 is quite lower than Fq. 

Fig. 3 Our results F2, i^s and F4 are compared with the the second and third-order results by the optimized 
expansions which were provided by Okopinska (z=0). In this figure, from the upper to lower (for large values of /3), 
the first, the second, the sixth and the seventh curves are the variational, the third-, the second- and the fourth-order 
results of our scheme, respectively. The third and the fourth curves are the second- and the third-order results of the 
optimized expansions respectively, and the fifth curve is the exact result. 

TABLE I. Our results F2, F3 and F4 are compared with the variational result Fq and the accurate free energies Faccu provided 
by Okopiriska (z=10). 



T 


F4 


-T accu 


Fo 


F2 


F3 


1. 


2.262259 


2.26225951564 


2.262452 


2.2622504 


2.262261 


2. 


2.063913 


2.06391575514 


2.064409 


2.0638734 


2.063925 


3. 


1.555676 


1.55569718863 


1.556991 


1.5555342 


1.555747 


4. 


0.7808495 


0.780936961496 


0.7836171 


0.7805129 


0.7811028 


5. 


-0.2099735 


-0.209722583045 


-0.2050294 


-0.210593 


-0.2093154 


10. 


-7.37775 


-7.37249823358 


-7.348171 


-7.3793287 


-7.367283 


20. 


-28.03925 


-27.9670036469 


-27.86147 


-28.0074342 


-27.92105 


30. 


-53.50769 


-53.2269143165 


-52.99767 


-53.3278138 


-53.08786 



TABLE IL Our results F2 and -fs are compared with the variational result Fq, the first- and third-order result, F^^ and 
F^ ^ obtained by the cumulant expansions ( F2 and F3 in Table II of Ref. [6]) and the exact free energies Fexa- In this table, 
m = LJ = 1 



A 


P 


Fo 


pKr 


Fs 


F.f^ 


Fexa 


F2 


1.0 


5.0 


0.812491 


0.81188 


0.807364 


0.803882 


0.803758 


0.800767 


5.0 


5.0 


1.244312 


1.24353 


1.2355 


1.22494 


1.22459 


1.216996 


50.0 


5.0 


2.54758 


2.54675 


2.529673 


2.50067 


2.49971 


2.480384 


500.0 


10.0 


5.425756 


5.42536 


5.387961 


5.32211 


5.3199 


5.276719 


20000.0 


3.0 


18.50166 


18.5003 


18.37314 


18.1449 


18.137 


17.98822 
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